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Detailed heat-transfer rate distributions measured laterally over the windward surface of an orbiter-like 
configuration using thin-film resistance heat-transfer gauges and globally using the newly developed relative inten- 
sity, two-color thermographic phosphor technique are presented for Mach 6 and 10 in air. The angle of attack was 
varied from 0 to 40 deg, and the freestream Reynolds number based on the model length was varied from 4 x 10 5 
to 6 X 10 A at Mach 6, corresponding to laminar, transitional, and turbulent boundary layers; the Reynolds number 
at Mach 10 w as 4 X 10 s , corresponding to laminar flow. The primary objective of the present study was to provide 
detailed benchmark heat-transfer data for the calibration of computational fluid -dynamics codes. Predictions from 
a Navier-Stokes solver referred to as the Langley aerothermodynamic upwind relaxation algorithm and an ap- 
proximate boundary-layer solving method known as the axisymmetric analog three-dimensional boundary layer 
code are compared with measurement. In general, predicted laminar heat-transfer rates are in good agreement 
with measurements. 


Nomenclature 

Cl, = heat-transfer coefficient, lbm/ft 2 -s, q/(h JW — where 

h\iw — hf 2 

C h ref - reference heat-translcr coefficient, lbm/ft 2 -s 

h = enthalpy, Btu/lbm 

L = model length in symmetry plane, in. (see Fig. 2) 

M - Mach number 

p - pressure, lb/in. 2 

q - heat-transfer rate, Btu/ft 2 -sec 

r = radius, in. 

Re - unit Reynolds number, ft 1 

T = temperature, R 

t = time, s 

V = velocity, ft/s 

,v = longitudinal distance, in. 
v = lateral or spanwise distance, in. 

a = angle of attack, deg 

y - ratio of specific heats 

Subscripts 

aw = adiabatic wall 

b = base 

L = based on model length in symmetry plane 

/, 1 = reservoir conditions 

/, 2 = stagnation conditions behind the normal shock 

li; = model surface; wall surface 

2 = static conditions behind the normal shock 

oo = freestream static conditions 

Introduction 

F OR many years, the Langley Research Center has been involved 
in the study of Earth-to-orbit space transportation concepts that 
fulfill a variety of anticipated mission needs. Some examples are pro- 
grams such as the Assured Crew Return Vehicle (ACRV), a vehicle 
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designed to return crew members from Space Station Freedom; the 
Personnel Launch System (PLS), a personnel carrier to low Earth or- 
bit and return; and Advanced Manned Launch Systems, a candidate 
replacement for the current Space Shuttle Orbiter (see tor example. 
Refs. 1-4). From an aerodynamic/aerothermodynamic perspective, 
these programs will benefit from knowledge gained as a result of 
the comprehensive data base established for the Space Shuttle pro- 
gram. Therefore, it is important that computational, ground-based, 
and flight data bases be brought together in an effort to gain an accu- 
rate knowledge of How field phenomena associated with the space 
transportation system. The impetus for this work is to provide im- 
provements to current computational fluid-dynamic techniques and 
to ground-to-flight extrapolation techniques that would be applied 
to the next space transportation system. 

The present study augments the well-established comprehensive 
aerothermodynamic data base for the Shuttle Orbiter by providing 
additional information concerning the complex three-dimensional 
windward flowfield for an orbiter-like configuration. Areas of inter- 
est include transition from laminar to turbulent boundary-layer 
heating phenomena and shock-shock interaction phenomena on 
windward-surface heating distributions. To date, most comparisons 
of aerothermodynamic results obtained via ground-based measure- 
ments, predictions, and orbiter flight data have focused on the wind- 
ward symmetry plane; similar comparisons for systematic spanwise 
heating are quite scarce. For the present study, detailed spanwise 
heat-transfer distributions were measured over the windward sur- 
face of a winged lifting entry configuration referred to herein as a 
modified orbiter. The objective of the present study is to compare 
predicted heating distributions from computational fluid-dynamics 
(CFD) codes with these detailed spanwise heating measurements. 
A set of high-fidelity modified orbiter models using the thin-film 
resistance gauge technique have been designed, fabricated, instru- 
mented, and tested in two wind tunnels of the Langley Hypersonic 
Facilities Complex to obtain detailed aerothermodynamic data over 
a wide range of test conditions. The discrete heat-transfer measure- 
ments obtained using thin-film gauges are augmented by global ther- 
mal mappings and qualitative heat-transfer mappings obtained with 
the recently developed relative-intensity, two-color thermographic 
phosphor technique. These measurement techniques nicely cap- 
ture the previously documented streak-heating phenomena resulting 
from the upstream effects of bow-shock-wing-shock interaction. 

Experimental Method 

Facilities 

The Langley 20-Inch Mach 6 Tunnel is a blowdown wind tunnel 
that uses dry air as the test gas. The air is heated to a maximum 
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temperature of 1088 R by an electrical resistance heater, and the 
maximum reservoir pressure is 525 psia. A fixed- geometry, two- 
dimensional, contoured nozzle with parallel side walls expands the 
How to Mach 6 at the 20-in. -square test section. This tunnel is 
equipped with a bottom-mounted model injection-retraction system 
capable of injecting heat-transfer models from a sheltered position 
to the nozzle centerline in less than 0.6 s. The run time for this 
facility varies from 2 to 10 min. A description of this facility and 
calibration results are presented in Ref. 5. 

The Langley 3 1 -Inch Mach 10 Tunnel, formerly known as the 
Langley Continuous-Flow' Hypersonic Tunnel, 6 is a blowdown faci- 
lity having a run time of approximately 60 s. The facility uses 
a water-cooled, three-dimensional contoured nozzle to generate a 
nominal Mach number of 10 at the 31 -in. -square test section. Dry 
air is used as the test gas and is heated to a maximum reservoir stag- 
nation temperature of 1900 R. The maximum reservoir pressure is 
1450 psia. The tunnel is equipped with a side- wall-mounted model 
injection-retraction system. Models are sheltered in the injection 
chamber on the side of the tunnel until the tunnel flow is started. 
After establishing Mach 10 flow at the test section, heat-transfer 
models are injected into the highly uniform stream within 0.5 s. 

Models 

Heat-transfer distributions were measured on the windward sur- 
face of a 0.0075-scale (9.09-in. length) model of a Shuttle-like 



a) Glass-ceramic model 



b) Stainless-steel model with glass-ceramic inserts 


Fig. 1 Photograph of modified orbiter. 



Fig. 2 Sketch of spanwise heating model. 



a) Complete geometry b) Modified geometry 

Fig. 3 Space Shuttle Orbiter geometry. Dashed lines represent fair- 
ings to model leeward side in present study. 

configuration. A photograph and a sketch of the model are presented 
in Figs. 1 and 2, respectively. The model has the same lower shape as 
the forward 93% of the Shuttle Orbiter. However, the upper surface 
is defined by elliptical cross sections. As shown in Fig. 3, the orbiter 
canopy has been faired and the vertical tail and OMS pods omitted. 

Two 9.09-in. -length models, referred to herein as the modified 
orbiters, were fabricated. Both models had solid steel upper sur- 
faces; however the lower surfaces of the two models differed. One 
model was machined entirely of steel and then slotted to accept 
machinable glass-ceramic substrates (herein referred to as the slot- 
ted stainless-steel model), whereas the second model was unique, 
being machined entirely from the ceramic material and slotted to 
accept the machinable glass-ceramic substrates (herein referred to 
as the slotted ceramic model). 

Both models were cut on a numerical milling machine using a 
tape generated with the geometry program described in Ref. 7. Since 
this program was also used to generate the geometry in the CFD 
codes (to be discussed subsequently), differences between the ex- 
perimental and computational models were within the machining 
tolerance of ±0.003 in. 

Test Conditions 

For this study, the pitot pressure could not be conveniently mea- 
sured when the model was positioned in the test-section of the 
31 -Inch Mach 10 Tunnel. Thus, test-section flow conditions were 
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Table 1 Nominal reservoir and freestream flow conditions 



Pt, 

psi 

T t . i 
R 

Po o» 

psi 

R 

10-5 

T w /T t , 2 

Oi, ret'* 

10- 5 Btu-s/ft 4 

5.82 

29 

872 

0.022 

1 12.2 

4.16 

0.62 

0.276 

5.94 

127 

915 

0.086 

1 13.6 

17.21 

0.59 

0.556 

5.98 

253 

911 

0.163 

111.7 

33.28 

0.59 

0.772 

6.01 

475 

931 

0.297 

1 13.2 

59.79 

0.58 

1 .050 

9.75 

354 

1822 

0.009 

94.9 

4.17 

0.30 

0.326 



based on measured reservoir pressures and temperatures and a recent 
unpublished calibration of the facility. 

Flow conditions in the 20-Inch Mach 6 Tunnel were determined 
from the measured reservoir pressure and temperature and the mea- 
sured pitot pressure at the test section. 

Nominal reservoir stagnation and corresponding freestream flow 
conditions for the present study are presented in Table 1. The value 
of L used to determine Re x L is 9.09 in. 

The stai nless-steel and ceramic models were tested over a range of 
angles of attack from 0 to 40 deg in 10-deg increments for each flow 
condition at Mach 6. Due to time constraints on testing in the 31- 
Inch Mach 10 Tunnel, only the stainless-steel model was tested at 
Mach 10; this model was tested at angles of attack equal to 0, 10, 15, 
20, and 30 deg. The sideslip angle was zero for all tests. The angle 
of attack was measured relative to the model centerline (see Fig. 2). 

Instrumentation and Testing Techniques 

77/ in - Fi l/n G auges 

Thin-fllm resistance heat-transfer gauges were used to mea- 
sure surface temperature-time histories from which heat-transfer 
rates were inferred. The technology of the thin-film gauges re- 
mains unchanged from that developed for the Langley Expansion 
Tube.* 9 The model surface contained 149 palladium gauges (0.04 
by 0.05 in.), each approximately 1000 A thick. These gauges were 
deposited on highly polished substrates that were precision-fitted 
to the model. The model substrate thickness was sized to provide a 
maximum run time of 1 .5 s; that is, the substrate essentially behaves 
as a semi-infinite slab for 1 .5 s over most of the forebody. An alu- 
minum oxide overlayer, approximately 5000 A thick, was deposited 
over the sensing elements as a means of increasing the gauge dura- 
bility. Detailed discussions of gauge construction, circuitry, and cal- 
ibration procedures and the data acquisition system are discussed in 
Refs. 8 and 9. The instrumentation layout is presented in Fig. 4. 

Thermographic Phosphor Technique 

The relative-intensity two-color thermographic phosphor tech- 
nique 111 12 is rapidly becoming the most widely used heat-transfer 


measurement technique in the Hypersonic Facilities Complex at 
Langley. 1 3 The model, generally M ACOR, stycast, or ceramic (fused 
silica), is coated with thermographic phosphors. When illuminated 
with ultraviolet light, electrons are excited and emit visible light dur- 
ing their subsequent relaxation to lower energy levels. The proba- 
bility that the relaxation and subsequent fluorescence emission, will 
occur is temperature-dependent, and the intensity of fluorescence 
may be used to determine local temperatures. Phosphor materials 
are selected for a two-color emission-band spectrum that matches 
front-end filters on conventional color video cameras. A three-chip, 
co-site sampling CCD camera, which has true color separation in the 
red, green, and blue bands and spatially congruent detection arrays 
is used. Early applications of this technique utilized thermographic 
phosphors having secondary and primary emission bands at 450 and 
520 nm, respectively (i.e., blue-green). Two materials were used, one 
with a temperature sensitivity range of 520 to 71 O R and the other 
with a range of 560 to 810 R. Currently, a mixture is used con- 
sisting of both a broadband phosphor and a narrowband rare-earth 
phosphor. For this mixture, the green and red camera filters are used, 
and a temperature sensitivity range of 480 to 860 R is obtained. The 
camera response is calibrated versus incident intensity. Measured in- 
tensities from two of the camera color outputs are used to form the 
intensity ratio from which quantitative temperature information may 
be determined using digital processing. Heating rates are calculated 
from surface temperature measurements using one-dimensional 
semi-infinite solid heat-conduction equations, as is discussed in de- 
tail in Refs. 1 1 and 12. Based on considerations presented in Refs. 12 
and 14, the uncertainty in the heating coefficients is ±15% . 

For this study, the heat-transfer model was cast from a high-purity 
fused-silica ceramic using an investment slip-casting technique de- 
scribed in Refs. 1 1 and 15. The slip casting forms a ceramic shell, 
which is then heat-treated to 2600 C R. A hydraulically setting mag- 
nesia ceramic is used to backfill the ceramic shell, thus providing 
strength and support to the sting structure. The phosphor mixture is 
suspended in a colloidal silica binder and spray-coated on the model 
surface. Typical coating thicknesses vary from 0.001 to 0.003 in. 

Data Reduction and Uncertainty 

The numerical method used to compute values of the heat-transfer 
rate q from the output of the thin-film resistance gauges is discussed 
in Refs. 9 and 16. A more accurate determination of the thermal 
properties of the glass-ceramic material 17 was used to reduce the 
data for the present study from those presented in Ref. 9. Sample 
time histories for gauges along the windward centerline revealed 
that the heat-transfer coefficient C h was essentially constant (i.e., 
within ±2%) over the time interval from 0.5 to 1.7 s; that is, the 
substrate behaved one-dimensionally. Second-order least-squares 
curve fits were applied to time histories of q for 0.5 < / < 1.7 
s, and values of q (or C h ) presented herein were obtained from 
these curve fits and generally correspond to / = l .1 s. For gauges 
along the planform leading edge, C h was observed to increase with 
time for t > 0.7 s, indicating that gauges in this region of small 
surface radius do not behave one-dimensionally. as expected. How- 
ever, for 0.5 s < / < 0.7 s, C h was essentially constant with time; 
thus values ol q (or C h ) for gauges at the planform leading edge cor- 
respond to r = 0.5 s. Heating distributions are presented in terms 
of the ratio of heat-transfer coefficients C h /C h , rc f. where C/, rct cor- 
responds to the stagnation-point heat-transfer rate to a sphere with 
a radius of 0.2067 in. For this 0.0075-scalc model, the nose radius 
(i.e., equivalent sphere radius) was determined to be 0.2067 in. This 
value was determined by scaling the approximate equivalent sphere 
radius for the full-scale Space Shuttle Orbiter nose, which is about 
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a) a = 15 deg 


c) a = 30 deg 
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Fig. 5 Global surface-temperature mappings of the windward surface of the modified orbiter using the two-color thermographic phosphor technique 

at Mach 10 and = 4.2 x 10 5 . 


27.6 in. (Ref. 1 8 ). Values of C hKf (see Table 1 ) were computed with 
the theory of Fay and Riddell 19 and correspond to a wall temperature 
of 540 R (that is, to an isothermal model). This method of presenting 
the measured heating (i.e., T w — 540 R) was possible because the 
measured heat-transfer coefficient was constant with time. 17 Run-to- 
run repeatability for the thin-film heat-transfer technique for oc = 0 
and a Reynolds number of 4. 16 x 10 s at Mach 6 was within ±1.1%; 
for Re^c L = 17.21 xlO 5 , a maximum run-to- run variation of ±3.5% 
was noted. Probable sources of error for thin-film resistance gauges 
arc discussed in Refs. 9 and 20. Based on these discussions and 
the present run-to-run variations, the heat-transfer rates presented 
herein are believed accurate to within ±8%> for all measured data. 

Prediction Methods 

AA3DBL Code 

Surface heating rates for the modified orbiter geometry were cal- 
culated using the axisymmetric analog for three-dimensional bound- 
ary layers developed by Cooke 21 and applied by Hamilton et al. 22 
The method is known as the AA3DBL, or axisymmetric analog 
three-dimensional boundary-layer, code. In using the AA3DBL 
code, several assumptions are made. As discussed in Ref. 22, the 
general, three-dimensional boundary-layer equations are written 


along a streamline. If the crossflow velocity in the boundary layer is 
neglected, the boundary-layer equations reduce to the axisymmet- 
ric form, provided that the distance along a streamline is interpreted 
as the distance along an equivalent body and the metric that de- 
scribes the spreading of the streamlines is interpreted as the radius 
of the equivalent axisymmetric body. This allows any axisymmet- 
ric boundary-layer solution to be used to calculate the approximate 
three-dimensional heating rates along streamlines. Thus, the heating 
over a complete vehicle can be computed by considering multiple 
streamline paths at a fraction of the cost that would be incurred if 
the full three-dimensional boundary-layer equations were used. 

Streamline information and boundary-layer-edge properties are 
acquired from the inviscid HALIS code. 23 No iteration is consi- 
dered between the boundary layer and inviscid solutions; thus, the 
Reynolds numbers must be reasonably large so that the boundary- 
layer solution does not strongly affect the inviscid solution. A rapid, 
approximate integral method 24 is used to calculate the heating along 
individual streamlines. For this study, HALIS/AA3DBL was exe- 
cuted for a perfect gas. The approximate run time for the HALIS 
code using a grid of 15 x 73 x 145 is 0.8 h on the Cray Y-MP 
computer for this configuration. The AA3DBL code requires ap- 
proximately 0.2 h on a Sun SPARC 2 using 145 streamlines. 2 * 5 
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LAURA Code 

Surface heating rales for the modified orbiter geometry were also 
calculated using the Langley aerolhermodynamic upwind relaxation 
algorithm ( LAURA). This code has been continually improved over 
the past few years and is described in detail in Refs. 26-29. (The 
latest version of the code is described in Ref. 29.) LAURA is a 
linite-volume-based algorithm that employs a point-implicit relax- 
ation procedure for obtaining the numerical solution to the govern- 
ing equations (Navicr-Slokes) for three-dimensional, viscous, hy- 
personic Hows, including chemical and thermal nonequilibrium. For 
the present study, LAURA was exercised for a perfect gas. The ap- 
proximate run time for the LAURA codeusing a gridof 70x81 x 145 
on the Cray-2 computer is 50 h for this configuration. With an in- 
crease in grid to 81 x 81 x 145, the approximate run time is 60 h. 

Results and Discussion 

Because of the large number of data obtained for the modified 
orbiter model in the two facilities, it is not possible to present all of 
the data in the limited space of this report. Thus, sample heating-rate 
distributions measured along the windward centerline and across the 
span arc presented herein. 

The modified orbiter was tested over a range of angles of at- 
tack at both Mach 6 and Mach 10 in air, and over a wide range of 
Reynolds numbers, corresponding to laminar, transitional, and tur- 
bulent heating-rate measurements, at Mach 6. Results were obtained 
at one value of Reynolds number at Mach 1 0 and correspond to lami- 
nar healing rates. The thermographic phosphor technique was used 
lo augment quantitative measurements at Mach 10. Data illustrating 
the effects of angle of attack and Reynolds number at Mach 6 were 
obtained with the slotted ceramic model. Data illustrating the effects 
of Mach number and angle of attack at Mach 10 were obtained with 
the slotted stainless-steel model, since this was the only model tested 
at Mach 10. At low Reynolds numbers and Mach 6, corresponding 
to a laminar boundary layer, differences between measured heating 
distributions using the two different models were negligible. 


Effect of Angle of Attack 

Global Surface-Temperature Mappings 

Global surface-temperature mappings of the windward surface 
of the modified orbiter at Mach 10 and Re x .L — 4.2 x 10 5 are 
presented in Fig. 5 for a = 15, 20, 30, and 40 deg. These mappings 
were obtained using the relative-intensity, two-color thermographic 
phosphor technique. The data presented correspond to the lime of 
arrival of the model at the nozzle centerline, which required ap- 
proximately 0.5 s. For consistency, the same color temperature scale 
50 F < T < 280 F was used for all angles of attack. (It should be 
noted that the temperature range for the thermographic phosphor 
used in this study is 75 F < T < 340 F.) 

As expected, the windward surface temperature increases with 
increasing a. Temperatures in the nose region and wing leading 
edge arc observed to exceed the upper limit of the chosen color- 
lempcraturc scale. However, for the regions of interest, namely, 
spanwise locations (i.c., wings), the selected temperature scale is 
quite adequate. 

Traversing the configuration in a spanwise direction from the 
centerline at approximately three-quarters of the model length for 
ot — 1 5 and 20 deg (Figs. 5a and 5b) reveals two distinct striations in 
temperature. The most inboard streak is the result of the downstream 
effects of the compression brought about by the discontinuous plan- 
form shape. As a result of this leading-edge compression, a stagnated 
region is produced, which spills over onto the windward surface (as 
a result of the highly complex three-dimensional fiowfield), causing 
a local enhancement in the heating over the windward surface. The 
outboard striation pattern is the result of the interaction of the bow 
shock with the wing shock. 

Further, density contours of the fiowlicld, obtained from an invis- 
cid LAURA solution (taken from Ref. 25) and reproduced in Fig. 6, 
illustrate the complex nature of the fiowfield in the vicinity of the 
wing. As observed in Fig. 6, the start of the wing shock occurs at 
x/ L ^ 0.7, and the downstream effect of this compression produces 
an abrupt increase in the local surface temperature (sec Fig. 5). 
Likewise, the interaction of the bow shock with the wing shock 



x/L 


Fig. 6 Density contours in plane near leading edge (taken from 
Ref. 25). 
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Fig. 7 Effect of angle of attack on windward centerline (y/L = 0.0) 
heating distributions. 


(x/L % 0.75 in Fig. 6) generates a similar increase in local surface 
temperature and hcncc local heating (to be discussed subsequently), 
as noted in Fig. 5. These locally enhanced heating phenomena have 
been observed from entry heat-transfer data for the orbiter wing 
lower surface 30 31 and from wind-tunnel tests of 0.0175- and 0.01- 
scale models of the orbiter at Mach 6, 8, and 10 using the phase- 
change paint technique. 32 ’ 33 As a is increased to 30 deg, the local 
enhancements in windward heating arc still disccrnable, although 
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not as distinct. With a further increase in a to 40 deg, these streaks 
appear to vanish. 

Longitudinal Heat-Transfer Distributions 

Measured heating distributions along the windward centerline of 
the modified orbiter using the thin-film gauge technique are pre- 
sented in Fig. 7 for Mach numbers of 6 and 10 and for a freestream 
Reynolds number based on length, Re^ of 4.2 x 10 5 . The varia- 
tion of heat-transfer coefficient ratio, C*/C*. re f, with wetted length 
x/L is presented for a range of angles of attack. At this Reynolds 
number, the flow along the windward symmetry plane is well be- 
haved and is indicative of a laminar boundary layer. As expected, 
the heat-transfer coefficient ratio increases with increasing angle of 
attack. In general, increasing a from 0 to 10 deg produces an in- 
crease in heating of approximately a factor of 3, whereas this same 
10-deg increase in a from 30 to 40 deg produces a much smaller 
increase of about 40% along the body. 

Spanwise Heat-Transfer Distributions 

Spanwise distributions of Ch/C h rci on the wing lower surface for 
a range of angles of attack are presented for x/L values of 0.75 and 
0.85 at Mach 6 and 10, in Figs. 8 and 9, respectively. These data 
correspond to Re^ L = 4.2 x 10 5 . For both Mach numbers and 
values of x/L y trends in spanwise heating vary appreciably with 
increasing a, indicating significant changes in the flowfield. 

Fora < 20deg, Q/C^ref increases immediately from the center- 
line in the direction of the wing leading edge ( v / L ~ 0.25) for both 




Fig. 8 Effect of angle of attack on spanwise distributions of heat- 
transfer coefficient ratio (solid symbols for « = 20 deg shown for clarity). 


Mach 6 and Mach 10. For these lower angles of attack, a two-order- 
of-magnitude increase in Q/CVret occurs between the centerline 
and the wing leading edge. One possible explanation for this pro- 
nounced increase in heat-transfer coefficient ratio in the spanwise 
direction is the inboard flow caused by the increase in pressure along 
the wing leading edge and planform discontinuity (i.e., discontinuity 
in planform shape caused by the junction of the double delta wing). 

For a > 20 deg, a nearly constant (i.e., more 2-D-like) heat- 
transfer coefficient ratio is observed for y/L < 0. 1 . For y/ L > 0. 1 , 
Ch/Ch'Kf increases in the direction of the wing leading edge, At 
these conditions, the flow along the symmetry plane exhibits par- 
allel streamlines (to be discussed subsequently); that is, increasing 
angle of attack serves to increase the pressure along the windward 
centerline, thus “straightening” the local flow. (It should be noted 
that according to the predictions of Weilmuenster, 23 the majority of 
the flowfield over the windward surface is supersonic for a < 30 
deg with the exception of a small subsonic pocket near the nose of 
the vehicle. However, for a > 30 deg, the flow is characterized by 
a large embedded subsonic region.) 

With increasing angle of attack, the difference in C h /C h re ( mea- 
sured at the centerline and at the wing leading edge decreases. For 
example, at or = 0 deg, there is approximately a two-order-of- 
magnitude increase in heating between the centerline and the wing 
leading edge (y/L ~ 0.25). However, at a = 40 deg, there is less 
than an order-of-magnitude increase in heating as measured at these 
two points. 
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b) x/L = 0.85 

Fig. 9 Effect of angle of attack on spanwise distributions of heat- 
transfer coefficient ratio (solid symbols for cx = 15 deg shown for clarity). 
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At Mach 10 (Fig. 9a), a local increase in heating is observed 
around y/L & 0.17 to 0.20 and x/L = 0.75 for a > 10 deg. At 
a — 30 deg, the local increase in C h /C h%K f is still evident but not 
as distinct as for a < 30 deg. (Thermal mapping is presented in 
Fig. 5c.) This discontinuity in spanwise heating trend is the result 
of increased pressure brought about by the discontinuous planform 
caused by the double delta wing and is discussed in the previous sec- 
tion entitled Global Surface Temperature Mappings; similar results 
are noted for x/L = 0.85 (Fig. 9b). 

Returning to the Mach 6 results of Fig. 8b, along with the first 
abrupt increase in C/,/C/, re f identified at y/L ~ 0.2 and a > 20 deg 
and discussed previously, a second well-defined discontinuity in the 
trend of the heat-transfer coefficient ratio occurs at y/L ~ 0.30 for 
a > 0 deg. For Mach 10 (Fig. 9b), this discontinuity occurs fora > 0 
deg and y/L 0.28. For both Mach numbers, asa is increased to 30 
deg or greater, this local increase in C;,/ Q >rc f along x/L = 0.85 is 
still visible, but not as pronounced. This region of locally enhanced 
heating is the result of the interaction of the bow shock with the wing 
shock. The use of closely spaced thin-film gauges allows sufficient 
spatial resolution to acquire discrete quantitative measurements in- 
dicating the magnitude of the streak-heating phenomena. 

Reynolds-Number Effects 

Longitudinal Heat-Transfer Distributions 

The effect of Reynolds number on centerline heat-transfer distri- 
butions at Mach 6 is presented in Fig. 10 for a range of Reynolds 
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Fig. 10 Effect of Reynolds number on longitudinal (y/L = 0.0) distri- 
butions of heat-transfer coefficient ratio at = 6. 


number from 4.2 x 1 0 5 to 59.8 x 10 5 and at angles of attack of 20 
and 40 deg. Differences between the measured data for the lami- 
nar cases [Re xL 5 33.3 x 10 5 for a = 20 deg (Fig. 10a) and 
Re xL S 1 7.2 x 1 0 5 for ot = 40 deg (Fig. 1 0b)] are essentially neg- 
ligible, (The effect of Reynolds number on laminar centerline heat- 
transfer coefficient ratio has been shown in Ref. 1 8 to decrease with 
increasing angle of attack for 1 5 deg < ot < 45 deg, corresponding to 
decreasing viscous effects.) Boundary-layer transition is observed 
to occur along the windward centerline around x/L r* 0.55 for 
Re x L = 59.8 x 10 5 and a = 20 deg and for Re x L = 33.3 x 10 5 
for a = 40 deg. For the highest Reynolds number and ot = 40 deg, 
the transition location has moved upstream toward the nose and is 
located around x/L ~ 0.25. The distribution of C*/CV rc f follow- 
ing this transition region is indicative of a fully turbulent boundary 
layer. At a =40 deg, the heating level for the highest -Re^x case 
(turbulent boundary layer) at x/L as 0.75 is about 5 times the lam- 
inar value. The value of turbulent heat-transfer coefficient at this 
location (x/L = 0.75) is very close to the value measured in the 
nose region (i.e., x/L = 0.075). 


Spanwise Heat-Transfer Distributions 

The effect of Reynolds number along rays normal to the center- 
line for a range of Reynolds numbers and a = 20 and 40 deg at three 
longitudinal locations (x/L = 0.35, 0.75, and 0.85) is presented in 
Figs. 1 1 and 12. For x/L = 0.35 (upstream of the planform dis- 
continuity) and ot = 20 deg (Fig. 1 la), C/,/Q let distributions are 
indicative of a laminar boundary layer for all values of Reynolds 
number, with essentially no effect of Reynolds number on spanwise 
heating. At a longitudinal station of 0.75 (i.e., downstream of the 
planform discontinuity) and a = 20 deg (Fig. 1 lb), heat-transfer 
coefficient distributions indicative of a laminar boundary layer are 
observed for Re 0 c L < 33.3 x 10 5 but not at the highest value of 
RtocX' Traversing across the span for Re xL = 59.8 x 10\ the 
value of Ch/ Ch n f at the centerline is indicative of a nonlaminar 
boundary layer followed by a return to laminar heat-transfer coef- 
ficient ratios prior to an abrupt discontinuity in the heating trend 
occurring at y/L ~ 0.20. The location (y/L & 0.20 for a = 20 
deg and Re Xt L = 59.8 x 10 6 ) of this discontinuity in Ch/Ch, n f cor- 
responds to the local downstream effects of the increased heating 
along the planform leading edge coinciding with the region of flow 
compression at the junction of the double delta wings. 

Continuing to the next longitudinal station at x/L — 0.85 and 
a = 20 deg (Fig. 1 lc), laminar heat-transfer coefficient ratios are 
observed adjacent to the windward centerline (i.e., y/L < 0.17) 
for Re x L < 33.3 x 10 5 , and evidence of local enhanced heating 
is observed for y/L locations of approximately 0.2 and 0.3. These 
abrupt increases in Ch/C hjtf have been cited previously as local 
downstream effects of increased heating along the wing leading 
edge, resulting from flow compression at the junction caused by 
the double delta wings (y/L 0.20) and the bow-shock-wing- 
shock interaction (y/L — 0.3). However, for Re x L > 33.3 x l(>\ 
the two distinct striations in heating may, in addition to the causes 
outlined above, be catalyst to local transition. At the wing leading 
edge and up onto the leeward surface, the effect of Reynolds number 
is essentially negligible for the laminar distributions. 

As Re oc i is increased above 33.3 x I0 5 , flow along the centerline 
(y/L = 0) is nonlaminar (probably transitional; see Fig. 10a) and 
then laminar again for 0.05 < y/L < 0. 15 away from the centerline, 
but then becomes transitional again at x/L ^ 0. 17. The flow again 
exhibits nonlaminar trends in C h /C h rc f , which are indicative of a 
fully turbulent boundary layer to the planform leading edge. These 
trends in spanwise heating were repeatable. The behavior of the flow 
along the span indicates that transition begins in regions or pockets. 

As discussed in Ref. 22, the manner in which transition occurs 
depends on the history of the streamlines producing the effect. For 
this reason, transitional and even turbulent flow could occur for one 
set of streamlines while the adjacent set of streamlines remained 
laminar. Qualitatively, this can be seen from the predicted inviscid 
surface streamline pattern shown for M x =9.15 and ot — 34.8 deg 
in the flight case in Fig. 1 3 (taken from Ref. 22). This flow behavior 
was also observed from other flight data 30,31 (Fig. 14). For a = 40 
deg and these same longitudinal locations (Fig. 12), similar heat- 
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Fig. 11 Effect of Reynolds number on span wise distributions of heat- 
transfer coefficient ratio: M ^ = 6; a = 20 deg. 


Fig. 12 Effect of Reynolds number on spanwise distributions of heat- 
transfer coefficient ratio at Mach 6 and a = 40 deg. 
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Kijr. 1 3 Computed surface streamline pattern on full-scale Shuttle Or- 
biter (taken from Kef. 22). 



Fig. 14 Transition front on Shuttle Orhiter wing from STS-2 flight 
(taken from Ref. 31). 
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Fig. 15 Effect of Mach number on longitudinal ( y/L = 0.0) distribu- 
tions of heat-transfer coefficient ratio. 


jin: trends are observed, although transition to a turbulent boundary 
layer is noted to occur at lower values of Re^x than for a = 20 deg. 
For this angle of attack and Re^x < 17.2 x I0\ the distributions 
of the heat-transfer coefficient ratio are well behaved and heating 
trends indicative of a laminar boundary layer are noted. As Re^x i s 
increased to 33.3 x 1()\ the C, f /CY re * distributions at x/L = 0.75 
and 0.85 become nonlaminar (transitional) at the centerline, with 
a general tendency to return to values characteristic of a laminar 
boundary layer in the spanwisc direction (i.c., increasing y/L). For 
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Fig. 16 Effect of Mach number on spanwi.se distribution of heat- 
transfer coefficient ratio at x/L = 0.85 and Re 0 0 x - 4.2 X 10 s . 

Re^x — 59.8 x 10\ the distributions imply a fully turbulent bound- 
ary layer with heat-transfer coefficient ratios approximately 4 to 5 
times their laminar counterparts. 


Effect of Mach Number 

Longitudinal Heat-Transfer Distributions 

The effect of Mach number on centerline nondimensional heat- 
transfer coefficient distributions at angles of attack of 10, 20, and 
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Fig. 17 Comparisons of windward centerline (y/L = 0.0) distributions 
of heat-transfer coefficient ratio obtained with thermographic phos- 
phors and thin-film gauge techniques. 

30 deg is presented in Fig. 15. These data correspond to Re^x = 
4.2 x 10 s and thus to laminar heating distributions. For the highest 
angle of attack investigated, the Mach 10 C/,/CV ret distributions are 
slightly higher than those for Mach 6 (approximately 5 to 8% near 
the trailing edge of the body). One possible contributor to this in- 
crease in heating with Mach number is the corresponding increase 
in boundary-layer displacement thickness, leading to a more blunt 
effective aeroline. Similar results are noted fora = 10 and 20 deg, 
with the difference between measurements at Mach 6 and 1 0 increas- 
ing slightly with decreasing or. Also, it should be noted that although 
the Reynolds number is matched closely between Mach 6 and 10, 
the ratio of wall temperature to total temperature is not; T w /T t 2 is 
approximately 0.3 at Mach 10 and 0.6 at Mach 6. The lower value 
of T w /T, 2 (i.e., colder wall) at Mach 10 may also contribute to the 
observed increase in heating with increasing Mach number. 

Spanwise Heat-Transfer Distribution 

The effect of Mach number on a spanwise distribution corres- 
ponding to a location x/L of 0.85 is presented in Fig. 16. Data 
are presented for a = 10, 20, and 30 deg at Re^ i — 4.2 x 10 5 . 
Again, on the centerline (y/L = 0.0) the Mach 10 data exceed 
the Mach 6 data by approximately 5 to 8%. This is also the case 
for 0 < y/L < 0.20. Beginning at y/L ^ 0.20 and continuing to 
the planform leading edge (y/L as 0.36), differences in location of 
enhanced heating produced by the leading-edge compression at the 
planform discontinuity and the local enhanced heating produced by 
the bow-shock-wing-shock interaction are observed. For example, 
at ot = 20 deg (Fig. 16b), where these llowfield characteristics are 
most prominent (see section entitled Global Surface-Temperature 
Mappings), the locations of these enhancements move farther in- 
board with increasing Mach number. At Mach 6, peaks in spanwise 
heating occur at approximately 0.22 and 0.30, whereas at Mach 
10 they occur at 0.20 and 0.26. This is expected and is a result 
of changes in shock structure (i.e., turning angles) with increasing 
Mach number. 

Comparisons with Two-Color Thermographic 
Phosphor Technique 

Ijtmxi tudinal Heat- Transfer Distributions 

Comparisons of surface heating distributions measured along the 
windward centerline using the thin-film gauge technique and the 
relative-intensity two-color thermographic phosphor technique are 
presented in Fig. 17. These data correspond to Mach 10, Re 0 Cii 
= 4.2 x 1 0\ and an angle of attack of 20 deg. For this laminar Dow 
condition, agreement between the two techniques is quite good, the 
maximum difference being less than 15%. The present comparisons 
are the first between the Langley standard for quantitative heat- 
transfer measurements (thin-film gauges) and the newly developed 
phosphor system. 10 ' 12 

Spanwise Heat- Transfer Distributions 

Comparisons of spanwnse heating distributions between the two 
measurement techniques are presented in Fig. 18 for a = 20 deg 
and x/L values of 0.65, 0.75, and 0.85. These data correspond to 


Re^.t = 4.2 x 10 5 . Again, the two techniques agree to within 
15%. It should be noted that the locations of the thin-film gauges 
are described in terms of a wetted- surface distance, whereas the 
thermographic phosphor technique uses a projected 2-D image of 
the windward planform. Thus, spatial differences exist between the 
two techniques, with the worst case occurring near the outboard 
wing section. At present the determination of exact longitudinal 
locations for spanwise comparisons is somewhat difficult and in- 
creases the inaccuracy, 10-12 The global information acquired using 
this technique is also an extremely important tool for deducing the 
complex 3-D flow-field effects on surface characteristics. Since mod- 
els cast for the thermographic phosphor technique arc orders of mag- 
nitude less expensive than those fabricated for use with the thin-film 
gauge technique and require significantly less time for fabrication, 
the thermographic phosphor technique is a very attractive means 
of obtaining near-quantitative global heat-transfer information in a 
rapid manner. 





Fig. 18 Comparisons of spanwise distributions of heat-transfer coeffi- 
cient ratio obtained with thermographic phosphors and thin-film gauge 
techniques. 
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Fig. 19 Comparison of measured and predicted heat-transfer rate dis- 
tributions along windward centerline ( y/L = 0.0). 





Fig. 20 Comparison of measured and predicted spanwise heat- 
transfer rate distributions. 


Comparisons with Prediction 

Longitudinal Heat-Transfer Distributions 

Sample comparisons of measured and predicted heating distribu- 
tions arc presented in this section for a Mach number approximately 
equal to 6 and Re x L ~ 4.2 x 1() 5 . A more complete set of theoreti- 
cal comparisons with the present experimental data set for angles 
of attack of 30 and 40 deg are presented in Ref. 25. Values of the 
heat-transfer rate q corresponding to a wall temperature of 540 R 
are used in these comparisons. 

Comparisons of measured and predicted heating distributions 
along the centerline are presented in Fig. 1 9 for an angle of attack of 
30 deg. For this laminar flow condition, the calculations are in very 
good agreement with measurement along the windward symmetry 
plane. Differences between the two computational methods and the 
experimental data are within approximately 10%. 

Spanwise Heat-Transfer Distributions 

Comparisons of measured and predicted heat-transfer rate to the 
wing lower surface are presented in Figs. 20a, 20b, and 20c for 
longitudinal locations x/L = 0.35, 0.75, and 0.85. respectively. 
These data are for a = 30 deg and correspond to laminar heat- 
transfer distributions. For x/L = 0.35 and 0.75, both the boundary- 
layer solver (AA3DBL) and the Navier-Stokes solver (LAURA) 
exhibit fairly good agreement with measurement over the wing and 
onto the leeward surface. For x/L = 0.85, LAURA provides a 
better qualitative description of the heat-transfer distribution, which 
includes the abrupt increase in heat-transfer rate at y/L & 0.20, 
than does the AA3DBL code. In general, both codes predict the 
spanwise heating levels reasonably well. 

Concluding Remarks 

Detailed distributions of the laminar, transitional, and turbulent 
heat-transfer coefficient ratio were measured laterally at several 
longitudinal stations over the windward surface of a modified Space 
Shuttle Orbiter configuration in air at Mach 6 and 10 for a range of 
Reynolds numbers and angles of attack. Predictions from an ap- 
proximate viscous method, AA3DBL, and a Navier-Stokes solver, 
LAURA, were compared with measured values in both longitudinal 
and spanwise directions. 

A significant increase in the measured heat-transfer coefficient 
ratio was observed with increasing angle of attack for both longitu- 
dinal and spanwise directions, as expected. For of < 10 deg, a two- 
order-of-magnitude increase was noted between the heat-transfer 
coefficient ratio measured at the centerline and in the region of 
the wing leading edge; this measured difference decreased with 
increasing angle of attack. With an increase in Reynolds number, 
the laminar flow over the windward surface became turbulent, and 
measured heat-transfer coefficient ratios for the turbulent flow cases 
were about 5 times those for the laminar flow cases. Spanwise dis- 
tributions indicate that transition occurs in regions or pockets; that 
is, for a particular longitudinal station the spanwise flow may be 
laminar and transitional combinations before becoming completely 
turbulent. Local enhancement regions in heating (i.e., striation or 
streak heating) over the wing were observed in spanwise distri- 
butions and are caused by wing leading-edge compression at the 
planform discontinuity producing a local stagnation region, which 
spills onto the wing lower surface, and by the interaction of the 
bow shock with the wing shock. The strength of each distinct stria- 
tion pattern decreased with increasing angle of attack. Heat-transfer 
rates predicted with the AA3DBL and LAURA codes were gener- 
ally in good agreement with measurement for both longitudinal and 
spanwise distributions. 
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